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TRACE FORMULAE AND INVERSE SPECTRAL THEORY 
FOR SCHRODINGER OPERATORS 



F. Gesztesy, H. Holden, B. Simon, and Z. Zhao 

Abstract. We extend the well-known trace formula for Hill's equation to general 
one-dimensional Schrodinger operators. The new function ^, which we introduce, is 
used to study absolutely continuous spectrum and inverse problems. 



In this note we will consider one-dimensional Schrodinger operators 

(IS) H = -— + V(x) on L^(R;dx) 

and Jacobi matrices 

(IJ) {hu){n) = u{n + 1) + u{n - I) + v{n)u{n) onf{Z). 

We will suppose that V{x) is continuous and bounded below and v(n) is bounded. 

In the analysis of the inverse problem for H when V is periodic (V{x + L) = 
V{x)), a crucial role is played by a trace formula [5, 13, 15]. H then has as 
its spectrum an infinite set of bands: spec(if) = [Eo,Ei] U [£'2, £^3] U Let 
{/i„(a;)}^]^ be the eigenvalues of the Dirichlet Schrodinger operator in L'^{x, x + L) 
(w.r.t. Lebesgue measure) with u{x) = u{x-\-L) = boundary conditions (i?2n-i < 
^^n{x) < E2n)- The trace formula says that if V is in _ff ^'^([0, L]), where is 
the Sobolev space of distributions with derivatives up to order m in L^, then 

oc 

(2) V{x)^Eo + J2{E2 n 

+ E2n-1 - 2fln{x)). 

n=\ 

One of our main goals here is to prove a version of this trace formula for arbitrary 
Schrodinger and Jacobi operators. 

We will need the paired half- line Dirichlet operator ifg defined on L^(— cxd, x) © 
L^(a;, 00) and /ip on /^(Z|to < n) ® P{Z\m > n) with ^(a;) (or u{n)) vanishing 
boundary conditions. In the periodic case, it can be shown that iJ,n{x) are precisely 
the eigenvalues of i?g (as long as -E2n-i < f^n{x) < £2™, i-e., no equality). 
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The difference {H — — (H^ — i)^^ is rank 1 (and similarly in the case of 
/ig if we define (/ig — i)~^{n,m) = 0)) and so trace class. As a result, the Krein 
spectral shift [11] exists; i.e., there is a function ^{x, A) uniquely determined a.e. in 
A w.r.t. Lebesgue measure by 

POO 

(3) Tr(/(iJ) - /(i/g)) = - / / '(A)^(x, A) dX, 



(4) 0<$(x,A)<l, 

A) = if A < inf (spec(iJ)) 

for any function, /, with sup;^ |(1 + A^) df /dX\ < oo. 

^ is a remarkable function which we claim is central to the proper understanding 
of inverse problems; it will be discussed in detail in three forthcoming papers which 
include detailed proofs of the theorems that we present here [6-8]. Our general 
trace formula is 

Theorem IS [6]. Let V be continuous at x and Eq < inf(spec(if)). Then 

/.oo 

(5S) V{x) = Eo + lim / e-"^(l - 2^{x, A)) dX. 

Theorem IJ [6]. Let E- < inf(spec(/i)) and E+ > sup(spec(/i)). Then 
(5J) v{n) = 1{E-+ E+) + j"^* Q - C(n, A)) dX. 

Remarks. 1. If F is smooth, there are higher-order trace relations including KdV 
invariants [7]. 

2. In the Jacobi case, ^(n. A) = 1 if A > sup(spec(/i)), which is needed for 
consistency in (5J). 

3. While we have singled out the Dirichlct boundary condition at x S R, any 
other selfadjoint boundary condition of the type il^^x) + /3tl){x) = 0, /3 G M, has 
been worked out as well in [7]. 

4. Besides the motivating equation (2), two other special cases are in the liter- 
ature. Kotani and Krishna [10] and Craig [3] discuss the case where V is bounded 
and continuous and (in our language) C = ^ a.e. on spec(iJ); and Venakides [16] 
has a trace formula when V is positive of compact support. In [6] we will discuss 
the relation of our work to these in more detail. 

Sketch of Proof. For simplicity, we consider only the Schrodinger case and suppose 
if > and take Eq = 0. By (3) 

/>OC 

Trie-"" - e-""^) = a er'^^^ix, X) dX. 
Jo 

Moreover, a path integral argument shows that 



TV(e-"^ - e-«^S) = i(i _ aV{x) + o{a)). 



Given that 



1 POO 1 

(6) 5=»/ e-'^- 

we get (5S) for Eq = 0. 

A second critical result that we prove is 
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Theorem 2 [6] . For each a; e R and a. e. A in M, 

^(x,A) = — &T:g{G{x,x]\ + iQ)). 

Remark. G is the integral kernel (resp. matrix elements) of (-ff — A)~^ (resp. [h — 
A)~^). By general principles for each x, lim^j^o G{x, x;X + ie) exists for a.e. A. 

Examples. 1. V = Q. In the H case, G{x,x;\) = (-A)-^/^ for A e C\[0,oo) 
with the branch of square root, so G > for A G (— oo,0). Thus, for A £ (0, oo), 
G{x,x;X + iO) = i|A|~^/^ and £,{x,\) = ^. Equation (6) is then an expression of 
the known fact that Tr{e-°'"° - e~""^-°) = \ for ah a. 

2. Let V be periodic and in iJ^'^([0, L]) with V{x + L) = V{x). The spectrum 
of H is U5^o[-^2n, -E2n+i] as noted already. Because V is in ff^'^([0, L]), 

oo 

(7) ^|£;2„--E2„-i| <oo. 

It can be shown (see, e.g., Kotani [9], Simon [14], and Deift and Simon [4]) that 
G{x, a;; A + iO) is pure imaginary on spcc(H), so ^ = ^ there. Thus wc claim (here 
and below, we do not give a value to ^ at points of discontinuity; the real-valued 
function ^ is only determined a.e.): 



$(x,A) 



E2n < A < -B2n+1) 

E2n+1 < A < Hn+\{x), 
Hn+lix) < A < E2n+2, 



for < 5 < 1, and ^ jumps by —1 at /i„+i(a;). Because of (7), \l — 2^{x, A)| rfA < 
oo and (5S) becomes (2). 

3. Let V{x) ^ oo as \x\ oo. Then H has eigenvalues Eq < Ei < E2 < ■ ■ ■ 
and iJg eigenvalues fii{x) < fi2{x) < ■ ■ ■ with i?„_i < /z„(a;) < En- |1 — 2^| = 1, so 
the integral in (5S) is not absolutely convergent if a is set equal to zero and (5S) 
becomes a summability result; explicitly 

00 

V(x) =Eo + lima-i X^pg-'^^ (^)« - e'^^" - e''^^-'"]. 

For an explicit case, let V{x) = x^ — 1 and place the Dirichlet condition at a; = 0. 
Then 

2n (n odd) 

2(n — 1) (n even, n > 2), 



En = 2n, /x„(0) 



so ^(0, A) = 1 on (0, 2) U (4, 6) U • • • and ^(0, A) = on (2, 4) U (6, 8) U • • • and 
formally 



/ 

Jo 



(1 -2^(0,A))dA = -2 + 2-2- 



'0 

The regularization (5S) is just the Abelian sum which is —1, which is exactly 1^(0). 
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4. Let V(x) he short range in the sense that V is L^(R). Then one can 
write down ^(x, A) in terms of the reflection coefficients -R(A) and Jost functions 
/+(x,A) (lim,^ooe-*^""V+(^,A) = 1), viz [8] 



1 1 



A > 0. 



l + i?(A)/+(a;,A)^ - 
l/+(^,A)|2 . 

In particular, |^(a;, A) - 5I < 5l-R(A)|, and if F e H'^''^{M.), we have that 



(8) ^(x,A) = - + -arg 

2 TT 



(9) / 

JEo 



dX < 00, 



so 



poo 

?o+ / (1- 



V{x)=Eo+ / (l-2C(a;,A))dA 



without a need for regularization. 

5. There is a general summability result [8] like (9) also for the sum of a smooth 
periodic potential and a sufficiently short-range potential modeling impurity scat- 
tering in one-dimensional crystals. 

The Krein spectral shift has rather strong continuity properties: 

Lemma 3a. Let Vm{x) {resp. Vm{n)) converge to V{x) uniformly for x € [—L,L] 
for each L {resp. to v{n) for each n) and so that 'mfx,mVm{x) < —00 {resp. 
suPn.m l^ni('^)l < ■ Tkcu as Tueasures in X, ^m{x,X)dX converges weakly to 

^{x, A) dX for each fixed x. 

It follows from Theorem 2 that 

Lemma 3b. For each fixed x, spec^^{H) = {A|0 < ^(A,a;) < 1}-°'^'' where is 
the essential closure. 

Third, it follows from results of Kotani [9] in the Schrodinger case and Simon 

[14] in the Jacobi case: 

Lemma 3c. IfV {resp. v) is periodic, then £,{x, X) = ^ onspec{H) {resp. spec(/i)). 

These three lemmas imply 

Theorem 3 [6]. Suppose Vm {resp. Vm) converge to V {resp. v) in the sense of 
Lemma 3a and each Vm {resp. Vm) is periodic. Then for any measurable set S CR 

15 n specie (i?) I > lim|S'nspec(ifm)| 
{resp. replacing H by h) where | • | = Lebesgue measure. 

Example. Consider the Jacobi matrix with v{n) = X cos(7ran) (almost Mathieu or 
Harper's model). Avron et al. [1] have proven that if a is rational, then | spec{ha)\ > 
4 — 2|A|. Theorem 3 then implies (by approximating any a by rationals) that 
\spec^^{ha)\ > 4 — 2|A|, slightly strengthening a recent result of Last [12]. In 
particular, we have a new proof of Last's spectacular result that spec^^^{ha) ^ if 
|A| < 2 and a is a Liouville number. 

Finally, [6] will use ^ to study the inverse problem. Typical of our results is the 
following: 

Let V{x) — > 00 as a; — *■ ±00. Let En{V) be the eigenvalues oi H = —d"^ /dx^ + V. 
We claim that when V is even, {En} are a complete set of spectral data in the 
sense that 
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Theorem 4. IfV,W are continuous functions onM. bounded from below, going to 
infinity at ±00, and obeying V{x) = V{—x) and W{x) = W{—x) so that En{V) = 
E„{W) for all n, then V = W. 

Borg [2] proved this result over forty years ago. The ^ function proof is natural, 
and we have an extension to the nonsymmetric case. When V is not symmetric, 
the Dirichlet eigenvalues and the information about whether each is a Dirichlet 
eigenvalue on (—00,0) or (0, co) also needs to be supplied. 
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